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Radio frequency waves do not penetrate into a plasma and are damped within it. The electric 
field of the wave and plasma current are concentrated near the plasma boundary in a skin layer. 
Electrons can transport the plasma current away from the skin layer due to their thermal motion. 
As a result, the width of the skin layer increases when electron temperature effects are taken into 
account. This phenomenon is called anomalous skin effect. The anomalous penetration of the rf 
electric field occurs not only for transversely propagating to the plasma boundary wave (inductively 
coupled plasmas) but also for the wave propagating along the plasma boundary (capacitively coupled 
plasmas) . Such anomalous penetration of the rf field modifies the structure of the capacitive sheath. 
Recent advances in the nonlinear, nonlocal theory of the capacitive sheath are reported. It is 
shown that separating the electric field profile into exponential and non-exponential parts yields an 
efficient qualitative and quantitative description of the anomalous skin effect in both inductively 
and capacitively coupled plasma. 



I. INTRODUCTION 



A radio frequency electromagnetic wave does not penetrate into a plasma if the wave frequency u> is smaller than 
the electron plasma frequency uj p = A-Ke 2 n e /m, where e and m are the electron charge and mass, respectively, and 
n e is the electron density. Electrons distribute their charge and current so as to shield out the electromagnetic wave. 
The shielding depends on the direction of the wave with regard to the plasma boundary. If the wave electric field is 
perpendicular to the plasma boundary, the rf field penetrates into the plasma only within a depth of the order of the 
Debye length vt/^ p , where vt — y/2T e /m is the electron thermal velocity, determined by the electron temperature 
T e , in eV. If the wave electric field is along the plasma boundary, the rf field penetrates into the plasma only within a 
depth of the order of the skin depth cjio v , where c is the speed of light in vacuum. Here, we consider a "collisionless" 
plasma, i.e. where the collision frequency is small compared to the wave frequency v <C u> and the electrons undergo 
rare collisions during the rf cycle; thus, collisions have little effect on wave screening by plasma. 

Another important scale is the nonlocality or phase-mixing scale vt/uj, which determines the scale length of the 
electron current profile in the plasma. To demonstrate the concept of phase-mixing scale Vt/<jJ let us consider a simple 
model, where an electron acquires a prescribed velocity kick near the plasma boundary, in the direction perpendicular 
to the boundary 

dv x (t) = AVexp(-iut). (1) 

The electron velocity at a distance x from the boundary will be determined by the moment when velocity kick was 
acquired at the plasma boundary, i.e., by the time t — x/v x . The electron current in the plasma is given by an 
integration over all electrons with a velocity distribution function f(v x ) 

/>OC 

j(x,t)=eAV f(v x )exp[-iu(t- x/v x )]dv x . (2) 
Jo 

Here, only electrons collided with the wall (v x > 0) have to be taken into account. For a Maxwellian distribution 
function f(v x ) = noe~ v */ VT /vt\/k, the plasma current in Eq.© becomes 

j{x,t) = ^— / exp -s H \ds, (3) 
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where s — v x /vx and jo — en AV. The amplitude and phase of the current are shown in Fig.l. In the limit 
lux/vt 3> 1, the integration in Eq. @ can be performed analytically making use of the method of steepest descend 
0, see Appendix A for more details. This gives 
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where X u — vt/ V%ui is the phase- mixing scale. Comparison of the asymptotic calculation result given by Eq. Q 
with the exact result of numerical integration in Eq. @ is shown in Fig. ^ From Fig. ^ it is evident that Eq. Q 
approximates the exact result for any x within a 15 percent error bar. The largest error occurs at x — 0, where half 
of the electron population with velocity v x > acquired the velocity kick, which gives rise to the electron current 
j(0) = jo/2, whereas Eq. (@J predicts j(0) — jo/V3, which corresponds to a 15 percent error. 




FIG. 1: Phase-mixing of the test particle current generated by velocity kicks AVcos(wi) at the plasma boundary: (a) current 
amplitude and (b) the current phase with respect to the phase of the velocity kick at the plasma boundary. The amplitude of 
the current is normalized on jo — enoAV, where no is plasma density. Solid lines show the exact result of numerical integration 
in Eq. dashed lines show the asymptotic, approximate analytical results given by Eq. @. 

Equation J3J describes the process of phase mixing - electrons with velocities different by Sv x ~ vt have different 
phase lag of the order ujx/vt at a distance x from the plasma boundary. Therefore, at x ~ vt /uj the phase difference 
becomes considerable: contributions to the total current from electrons with different velocities v x cancel out each 
other, and the plasma current vanishes. Interestingly, the spatial profile of the current is not a simple exponential 
function, but an exponential function of (x/X^) 2 ^ 3 . As it will be shown below this is typical for the spatial profiles of 
the electric field and electron current in warm plasmas due to nonlocal effects. 

So far, we solved only test-particle problem and did not take into account the plasma polarization. The current in 
Eq. J2J is nonuniform; thus, there must be an electron density perturbation according to the continuity equation 

dn e dj 

e ^ = -Tx (5) 

The electron density perturbations polarize the plasma and generate an electric field, which in turn, affects the electron 
motion and the electron current profile. Thus, Eq.© has to be modified to include the self-consistent electric field. 
This requires solving the Vlasov equation together with the Poisson equation. In his famous 1946 paper, Landau 
obtained an analytic solution for the penetration of the longitudinal rf electric field into a plasma |2j . Note that he 
also described "Landau damping" in the same paper. We briefly review his solution for a small amplitude electric 
field in the linear approximation and discuss the more realistic case of a large amplitude electric field. 

The structure of this review is as follows: In section II, the penetration of the longitudinal electric field into the 
plasma is described. This case corresponds to a capacitively coupled plasma. In section III, the penetration of the 
transverse electric field into the plasma is studied, which corresponds to an inductively coupled plasma. In subsection 
III.E, it is shown that anisotropy of the electron velocity distribution function can have a profound effect on the 
anomalous skin effect. 
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II. PENETRATION OF THE RF ELECTRIC FIELD DIRECTED PERPENDICULAR TO THE 
PLASMA BOUNDARY (CAPACITIVELY-COUPLED PLASMA) 

A. Small-amplitude electric field 



In the previous section, we considered a test particle current driven by artificially applied velocity modulations at the 
plasma boundary. Here, self-consistent penetration of a small amplitude rf electric field directed perpendicular to the 
plasma boundary is considered. Such a model provides some insight into the sheath structure of capacitively-coupled 
plasmas. 

First, let's consider a stationary negatively biased electrode. It is well-known that the externally applied el ectric 
field penetrates inside the plasma over distances of the order of the Debye length a — vt/ V2oj p = ^Te/^ne 2 ^. The 
plasma electrons are trapped by the plasma potential, <p(x), in the potential well ~e<j){x). The electron density obeys 
the Boltzmann distribution 



n e (x) = n exp [e<p(x)/T e 



(6) 



The Poisson equation 



d 2 ^ _ 
dx 2 

can be simplified assuming small potential variations 
fii = no- Thus, Eq. Q becomes 



-47re(n.i - n e ) (7) 
■e<f>{x)/T e <C 1 and a uniform background plasma with n e = 



d 2 (f> 47re 2 n 
dx 2 = T e 

The solution of Eq. (JSJ is an exponentially decaying electric field E — —d<f>/dx 

E = E exp(-- > ) . 



(8) 



(9) 



Here, Eq is the value of the electric field at the plasma boundary. This is the solution for a steady state, time- 
independent sheath electric field. In the opposite case of the time-dependent electric field, the Boltzmann distribution 
given by Eq. © is no longer valid and the electron density has to be determined from the Vlasov equation. Landau 
solved the Vlasov equation coupled with the Poisson equation analytically in the linear approximation considering an 
electrostatic wave with small amplitude \e(j>(x)\/T e -C 1 and small frequency u> <C co p |2|. Details of the solution are 
described in Appendix B. 

To summarize, the solution can be separated into three parts, 



E x {x,t)= £ cxp +E b + E t (x) 



(10) 



Here, Eq is the amplitude of the electric field at the plasma boundary, Eb = E^/e is the electric field in the plasma 
bulk far away from the sheath region, e = 1 — ui 2 /ui 2 is the dielectric constant of the cold plasma, and E t (x) is the 
electric field in a transient region with a spatial length of order vt/oj- The first term is the Debye screening of the 
external electric field. The second part describes a small, uniform electric field penetrating into the plasma far away 
from the boundary. The second and third terms are absent for a stationary applied electric field and appear only in 
the case of the rf electric field. The solution for the transient electric field E t (x) profile is derived in Appendix B and 
is given by 



E t (x) = / -- 

T JO k 1 



Im[eu(u}, k)] 
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where £\\(ui, k) is the longitudinal plasma permittivity (E||k) 
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and Z{Q is the plasma dispersion function 



Z(0 = tt- 1/2 / dt 



_ exp (— t 2 



-, 7mC > 0. 
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FIG. 2: Penetration of the external electric field into a plasma. Only part of the electric field E t (x) is shown. Solid lines 
show the exact solution given by Eq. Ijlljl : dashed lines and subscript appr correspond to the approximate calculation of Eq. ( 
1141 : dotted lines and the subscript std correspond to the approximate calculation in the limit x S> A„ obtained making use 
the method of the steepest descend given by Eq. p5^ . The rf electric field frequency is 13.56 MHz and the plasma density is 
10 8 cm" 3 (lines) and 10 9 cm" 3 (symbols). 



In the limit x 3> A w only small k contribute to the integral and £||(w, k) can be substituted by ejj(a>,0) = e in the 
denominator of Eq. I|14[) . which gives 

E t {x) « E appr (x) = — £ / -Imle^^^dk, (14) 
ire Jq k 

Application of the method of steepest descend to Eq. i|14|) yields 



, ,2 / \ 2/3 
2E w o / x ^ 



EAx) « E std (x) = —= §r — exp 
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where A w = vt/V^uj is the phase-mixing scale. The plots of the amplitude and phase of the electric field profile Et(x) 
given by Ea. l|ll|) and the approximate E appr (x) given by Eq. I|14l) . and asymptotic analytical result E st( i{x) given by 
Eq. H15|) are shown in Fig. El Figure |21 shows that the steepest descend method given by Eq. I|15l) closely approximates 
Eq. (|14|) already for x > vt/w. However, the both asymptotic solutions in Eq. (|14|) and Ea. H15|) approximate the full 
solution in Eq. (|llfl only for very large x > 40Vt/u- This is due to the made substitution e(tu,k) by e(w,0), which 
results in a considerable error for k ~ uj/vt or x ~ vt/oj. 

It follows from Eq. U5\) that the electric field amplitude at x > vt/w is of order Eq/s, i.e., it is comparable with 
the electric field far away from the boundary (E t ~ Ef, ). 

The origin of the electric field Et{x) can be explained by analyzing the individual electron dynamics. After passing 
through the region of the rf field, an electron acquires changes Ae(v x ) in energy and Au(v x ) in velocity 

/oo 
v x eE[x(t),t]dt, (16) 
-OO 

Au(v x ) = . 

mv x 

Here, the electron trajectory is x(t) = v x t, v x — \v x \sgn(t), and the electric field profile is given by Eq. lfTU|) . The 
total velocity kick is the summation over velocity kicks due to exponential, bulk and transitional electric fields 

Au(v x ) = Au + Au b + Au t . (17) 

Substituting an exponential electric field into Ea. (|16f) gives the corresponding electron velocity kick 

A , . 2eE Lo 2 a 2 , 

Au (v x ~ " . 18 

muj v x + [Loay 
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Substituting the uniform electric field Et, into Eq. (|16|l . gives the electron velocity 



Au b (v x 



2eE 



imcue 



(19) 



This calculation can also be explained as follows: An electron has the oscillating velocity Au s — eE^i/rruv in a uniform 
rf electric field and a thermal velocity v x . After a collision with the wall, an electron changes its velocity direction. 
If the initial average velocity was v x < 0, after the collision with the wall with specular reflection, the new average 
velocity v' x > will change according to 



or the average velocity changes to 



•y x + Au s (t) = -[v x + Au s (t)] 



-v x - 2Au s (t), 



(20) 



(21) 



which results in the effective velocity kick of Eq. 11911 . 

The origin of the electric field in the transition region E t {x) is due to the plasma polarization. The velocity 
perturbations Au s (v x ,t) produce bunches in the electron density, which, in turn, generate the electric field E t (x). 
The decay of the electric field E t (x) is due to phase mixing similarly to the test-particle case in Eq. Thus, 
generation of the transitional electric field E t (x) can be considered as a plasma self-consistency effect. 

The electric field E t (x) generates a significant portion of the total velocity kick and thus noticeably influences the 
electron heating in the rf electric field. Figure |21 shows the amplitude of the electron velocity kick Au(v x ) due to the 
interaction with the electric field given by Eq. (|10fl . Electrons with small velocities v x ~ u>a — vtuj /w p pick up a large 
velocity kick due to the exponential electric field Eq cxp(— x/a — iujt), Au ~ Auo ~ 2eEo/mu>. For very large electron 
velocities v x 3> vt, the velocity kick given by Eq. (|18fl becomes small and the main contribution to the velocity 
kick comes from the uniform electric field Et, = E^er wt je and the collision with the wall, Au ~ Aut, ~ 2eEo/moj£. 
In the intermediate range of velocities v x ~ vt, the account of the electric field E t (x) is important, as in this case 
Au t ~ Aut,. As is evident from Fig. El taking this electric field E t {x) into account results in a considerable reduction 
of the electron velocity kick v x ~ vt for the bulk of the electron population compared with the case when this electric 
field is not taken into account. Note that most models neglect the electric field E t (x), see for example 



LU 

a 
E 
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FIG. 3: Electron velocity kick after interaction with the rf electric field. Solid line shows a velocity kick Au calculated according 
to the full electric field in Eq. HUH . The dashed line shows a Am due to the electric field -Bo exp(— x/a) and the uniform electric 
field Eb only; the dotted line is due to Et- The rf electric field frequency is 13.56 MHz and the plasma density is 10 8 cm -3 . 



B. Large amplitude electric field 



In many practical applications, the value of the external electric field is large: the potential drop in the sheath 
region V s h is typically of the order of hundreds of Volts and is much larger than the electron temperature T e , which 
is of the order of a few Volts; consequently the electric field penetration has to be treated nonlinearly. 



6 



In the limit V s h T e , a wall is charged negatively all time with an alternating charge in a manner to conduct an 
ac current, driven by an external electric circuit. A negative charge pushes electrons away from the electrode up to a 
distance where the negative electric field is screened by a positive ion density. As V s h ^> T e , the sheath width is much 
larger than the Debye length and the plasma sheath boundary can be considered as infinitely thin. The position of 
the boundary is determined by the condition that the external electric field is screened in the sheath regions when 
and where electrons are absent (fj 0] • 

Electron interactions with the sheath electric field are traditionally treated as collisions with a moving potential 
barrier (wall). It is well known that multiple electron collisions with an oscillating wall result in electron heating, 
provided there is sufficient phase-space randomization in the plasma bulk. It is common to describe the sheath 
heating by considering electrons as test particles, and neglecting the plasma electric field @- As was pointed out in 
Refs. H I i accounting for the electric field in the plasma reduces the electron sheath heating, and the electron 
sheath heating vanishes completely in the limit of uniform plasma density. Therefore, an accurate description of the 
rf fields in the bulk of the plasma is necessary for calculating the sheath heating. The electron velocity is oscillatory in 
the sheath, and as a result of these velocity modulations, the electron density bunches appear in the region adjacent 
to the sheath, similar to the previously described case of small-amplitude wave, see Fig. These electron density 
perturbations decay due to phase mixing over a length of order Ux/w, where vt is the electron thermal velocity, and 
u> is the frequency of the electric field. The electron density perturbations polarize the plasma and produce an electric 
field in the plasma bulk. This electric field, in turn, changes the velocity modulations and correspondingly influences 
the electron density perturbations. Therefore, electron sheath heating has to be studied in a self-consistent nonlocal 
manner assuming a finite-temperature plasma. 




X 



FIG. 4: Schematic of density bunches formation in the region adjacent to the sheath. The plasma-sheath boundary is shown 
by bold solid line. Electrons with the same velocity v x and distance apart ~ v x /u collide with the sheath. The first electron 
looses its energy and decelerates, whereas the second acquire energy and accelerates. As a result, the distance between two 
electrons decreases, which produces electron density perturbations. 



Notwithstanding the fact that particle-in-ccll simulations results have been widely available for the past decade 
[Tot fTl| . a basic understanding of the electron heating by the sheath electric field is being incomplete, because no one 
has studied the electric field in the plasma bulk using a kinetic approach, similar to the anomalous skin effect for the 
inductive electric field @- In this regard, analytical models are of great importance because they shed light on the 
most complicated features of collisionless electron interactions with the sheath. In Ref.^|j an analytical model was 
developed to explore the effects associated with the self-consistent non-local nature of this phenomenon. 

One of the approaches to study electron sheath heating is based on a fluid description of the electron dynamics. 
For the collisionless case, closure assumptions for the viscosity and heat fluxes are necessary. In most cases, the 
closure assumptions are made empirically or phenomenologically [TT| . |13|. The closure assumptions have to be 
justified by direct comparison with the results of kinetic calculations as is done, for example, in Refs. [Til Il5| . 
Otherwise, inaccurate closure assumptions may lead to misleading results as discussed below. To model the sheath- 
plasma interaction analytically, the following simplifying assumptions have been adopted in Ref. [T^| . The discharge 
frequency is assumed to be small compared with the electron plasma frequency. Therefore, most of the external 
electric field is screened in the sheath region by an ion space charge. The ion response time is typically larger than 
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FIG. 5: Schematic of a sheath. The negatively charged electrode pushes electrons away at different distances depending on the 
strength of the electric field at the electrode. Shown are the density and potential profiles at two different times. The solid line 
is at the time of maximum sheath expansion. 



the inverse discharge frequency, and the ion density profile is quasi-stationary. There is an ion flow from the plasma 
bulk towards the electrodes. In the sheath region, ions are being accelerated towards the electrode by the large sheath 
electric field, and the ion density in the sheath region is small compared with the bulk ion density. In the present 
analytical treatment, the ion density profile is assumed fixed and is modelled in a two-step approximation: the ion 
density rib is uniform in the plasma bulk, and the ion density in the sheath n s h < rib is also uniform (see Fig.EJ. At 
the sheath-plasma boundary, there is a stationary potential barrier for the electrons (e$ s h), so that only the energetic 
electrons reach the sheath region. The potential barrier is determined by the quasineutrality condition, i.e., when the 
energetic electrons enter the sheath region, their instantaneous density is equal to the ion density [n e ((f> s /j) = n s h\. 

The electron density profile is time-dependent in response to the time- varying sheath electric field. The large sheath 
electric field does not penetrate into the plasma bulk. Therefore, the quasineutrality condition holds in the plasma 
bulk, i.e., the electron density is equal to ion density, n e — rib- In the sheath region, the electrons are reflected by the 
large sheath electric field. Therefore, n e — n s h for x > x s h{t), and n e — for x < x s h(t), where x 8 h{t) is the position 
of the plasma-sheath boundary [jj. From Maxwell's equations it follows that V • J = 0, where the total current J is 
the sum of the displacement current and the electron current. In the one-dimensional case, the condition V • J = 
yields the conservation of the total current 0, Q : 

1 dE h 

en e V e + - -f- = josm(u)t + <f>), (22) 

47r at 

where jo is the amplitude of the rf current controlled by an external circuit and <fi is the initial phase. In the sheath, 
electrons are absent in the region of large electric field, and Eq. (1221) can be integrated to give Q 

E sh (x,t) = -^[-1- cos(uit + 4>)] + 4Tr\e\n sh x, x<x sh {t) (23) 
u> 

where Poisson's equation has been used to determined the spatial dependence of the sheath electric field. The first 
term on the right-hand side of Eq. I|23|) describes the electric field at the electrode and the second term relates to 
the ion space charge screening of the sheath electric field. The position of the plasma-sheath boundary x s h(t) is 
determined by the zero of the sheath electric field, E g^x s h{t) , t) = 0. From Eq. (|23[l it follows that 

x sh (t) = — [1 + cos{ujt + (f)] , (24) 



where V s ho — jo/( en sh) is the amplitude of the plasma-sheath boundary velocity. The ion flux on the electrode 
is small compared with the electron thermal flux. Because electrons attach to the electrode, the electrode surface 
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charges negatively, so that in a steady-state discharge, the electric field at the electrode is always negative, preventing 
an electron flux on the electrode. However, for a very short time ((Jt n + w tt(1 + 2n)) the sheath electric field 
vanishes, allowing electrons to flow to the electrode for compensation of the ion flux. Note that there is a large 
difference between the sheath structure in the discharge and the sheath for obliquely incident waves interacting with 
a plasma slab without any bounding walls. Because electrodes are absent, electrons can move outside the plasma, 
and the electric field in the vacuum region, E s h(x,t) — (Airjo/ui) cos(wi + 0), may have an alternating sign. Therefore, 
electrons may penetrate into the region of large electric field during the time when E s h(x,t) > [la. Il7|. In the 
discharge, however, because the sheath electric field given by Eq. ll'.'iH always reflects electrons, the electrons never 
enter the region of the large sheath electric field, which is opposite to the case of obliquely incident waves. 

The calculations based on the two-step ion density profile model are known to yield discharge characteristics in 
good agreement with experimental data and full-scale simulations [Isj . 

For analytical calculation of the rf electric field inside the plasma, a linear approximation is used for the plasma 
conductivity. The validity of the linear approximation is based on the fact that the plasma-sheath boundary velocity 
and the mean electron flow velocity are small compared with the electron thermal velocity, V s h <C «r, 0- U3 The 
important spatial scale is the length scale for phase mixing, A w . The sheath width satisfies 2V s m/u -C A w because 
V s h *C Vt- Therefore, the sheath width is neglected, and electron interactions with the sheath electric field are treated 
as a boundary condition. The collision frequency (y) is assumed to be small compared with the discharge frequency 
(y -C u>), and correspondingly the mean free path is much larger than the length scale for phase mixing. Therefore, 
the electron dynamics is assumed to be collisionless. The discharge gap is considered to be sufficiently large compared 
with the electron mean free path, so that the influence of the opposite sheath is neglected. The effects of a finite gap 
width have been discussed in Refs. [T^.l20l|. 

The electron interaction with the large electric field in the sheath is modelled as a collision with a moving oscillating 
rigid barrier with velocity V s h{t) — dx s h{t) / dt. After a collision with the plasma-sheath boundary - modelled as a 
rigid barrier moving with velocity V s h(t) - an electron with initial velocity — u acquires a velocityit + 2V s h- Therefore, 
the power deposition density transfer from the oscillating plasma-sheath boundary is given by p 



Psh = j{J v du i u+V ^(-t)} [(2V ah (t)+u) 2 -u 2 ] f sh (- U) t)j, (25) 

where m is the electron mass, f s h(—u, t) is the electron velocity distribution function in the sheath, and (■ • •) denotes a 
time average over the discharge period. Introducing a new velocity distribution function g(—u', t) = f s h[—u — V s h{t), t], 
Eq. US) yields 



P sh = -2m(v sh (t) J u' 2 g{~u',t)du'^ , (26) 

where —u' = —u — V s h is the electron velocity relative to the oscillating rigid barrier. From Eq. l|2()|) it follows that, 
if the function g(u') is stationary, then (P s h — 0) and there is no collisionless power deposition due to electron 
interaction with the sheath |5l Il3t l21|. For example, in the limit of a uniform ion density profile n s h — rib, s(w') is 
stationary (in an oscillating reference frame of the plasma-sheath boundary), and the electron heating vanishes 
Q. Indeed, in the plasma bulk, the displacement current is small compared with the electron current, and from 
Eq. IL'l't it follows that the electron mean flow velocity in the plasma bulk, Vb(t) = —jo sm(cut + 4>)/enb, is equal to 
the plasma-sheath velocity V s h{t), from Eq. Ij24)). Therefore, the electron motion in the plasma is strongly correlated 
with the plasma-sheath boundary motion. From the electron momentum equation it follows that there is an electric 
field, Eb = m/edVb(t)/dt, in the plasma bulk. In a frame of reference moving with the electron mean flow velocity, 
the sheath barrier is stationary, and there is no force acting on the electrons, because the electric field is compensated 
by the inertial force (eEb — mdVb(t)/dt = 0). Therefore, electron interaction with the sheath electric field is totally 
compensated by the influence of the bulk electric field, and the collisionless heating vanishes Q . The example of a 
uniform density profile shows the importance of a self-consistent treatment of the collisionless heating in the plasma. 
If the function g(u',t) is nonstationary, there is net power deposition. In Ref. [T^. a kinetic calculation is performed 
to yield the correct electron velocity distribution function g(u',t) and, correspondingly, the net power deposition. 

The electron motion is different for low-energy electrons with an initial velocity in the plasma bulk |u| < u s h, where 
u 2 sh = 2e$> s h/m and for energetic electrons with velocity |m| > u s h. The low energy electrons with initial velocity — u 
in the plasma bulk are reflected from the stationary potential barrier e& s h, and then return to the plasma bulk with 
velocity u. High energy electrons enter the sheath region with velocity u\ = — (u 2 — u 2 ^) 1 / 2 . They acquire a velocity 
U2 = 2V s h — u\ after collision with the moving rigid barrier, and then return to the plasma bulk with a velocity 

(u 2 +u 2 sh )i/ 2 m- 
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As the electron velocity is modulated in time during reflections from the plasma-sheath boundary, so is the energetic 
electron density (by continuity of the electron flux). This phenomenon is identical to the mechanism of klystron 
operation pflj . The perturbations in the energetic electron density yield an electric field in the transition region 
adjusted to the sheath, see Fig0] 

The solution for the electric field E t (x) was obtained analytically in Ftef. [12j| . Similar to the previous section, the 
solution is an expression for the inverse Fourier transform. It cannot be represented in an analytical form and has to 
be simulated numerically. This simulation has been performed for n s h/rib — 1/3, lo/u> p = 1/100, and a Maxwellian 
electron distribution function. The electric field profile is close to E t (x) « E t o exp(— x/X c ), where E t o = — 0.72T e /A w , 
and A c = (0.19 + 0.77z)A w for x < 6Vr/w. For x > 6Vr/w, the electric field profile is no longer a simple exponential 
function, which is similar to the case considered in the previous section. The difference in phase of the currents of 
the energetic and low-energy electrons was observed in Refill]], but it was misinterpreted as the generation of the 
electron acoustic waves. Electron acoustic waves can be excited if there is a complex value of k, with small damping 
Im(k) -c Re(k), which is the root of the plasma dielectric function e(oj,k) = for a given ui. For a Maxwellian 
electron distribution function, such root does not exist when u> <C u> p . However, the electron acoustic waves can exist 
if the plasma contains two groups of electrons which have very different temperatures j^. The wave phase velocity 
is uj/k = y/ nc/nhy/Th/m , where n c and rih are the electron densities of cold and hot electrons, respectively, and 
T h is the temperature of the hot electrons. The electron acoustic waves are strongly damped by the hot electrons, 
unless n c <C and T c -C Th , where T c is the electron temperature of the cold electrons |24| . In the opposite limit, 
n c > 4rih, the electron acoustic waves do not exist pij . In capacitively-coupled discharges, the electron population 
does strati fy in to two populations of cold and hot electrons, as has been observed in experiments [23 and simulation 
studies [2a, |2y| . Cold electrons trapped by the plasma potential in the discharge center do not interact with the large 
electric fields in the sheath region and have low temperature. Moreover, because of the nonlinear evolution of plasma 
profiles, the cold electron density is much larger than the hot electron density [2{|. Therefore, weakly-damped electron 
acoustic waves do not exist in the plasma of capacitively-coupled discharges. Reference jll| used the fluid equation 
and neglected the effect of collisionless dissipation, thus arriving at the incorrect conclusion about the existence of 
weakly-damped electron acoustic waves. 




FIG. 6: Plot of the averaged square of the dimensionless velocity kick as a function of the dimensionless velocity for the 
conditions in Fig.l, taking into account (a) both E\{x) and Eb - solid line; (b) only Eb - dashed line; and (c) no electric field - 
dotted line. 

The power deposition is given by the sum of the power transferred to the electrons by the oscillating rigid barrier 
in the sheath region and by the electric field in the transition region, 

Ptot = Psh + Ptr- (27) 
Note that Ptr can be negative. Calculations making use of the Vlasov equation yield (l2^ 

Ptot = — [ muD u {u)'^-du, (28) 
,/n du 



where 



_ . „ u\du\ 2 

D u (u) = ^-^ (29) 
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is the diffusion coefficient in velocity space, and du is the change in the electron velocity after passing through the 
transition and sheath regions, 



du = 2iVb 



u n b 

0(u 

u n sh 



Ush) 



1 



eE t (k — lu/u) 



(30) 



where E t (k) is the Fourier transform of the electric field E t {x). First term describes the velocity acquired by fast 
electrons (|u| > u s h) in collisions with the sheath; the second is due to the bulk electric field Et, and collisions with 
either the potential barrier & s h or sheath; and the third is due the electric field in the transitional region E t (x). A 
plot of |du| 2 /4 is shown in Fig.|SJ Taking into account the electric field in the plasma (both Eb and E t ) reduces \du\ 
for energetic electrons (u > u s h) and increases \du\ for slow electrons (it < u s h)- Therefore, the electric field in the 
plasma cools the energetic electrons and heats the low- energy electrons, respectively. Similar observations were made 
in numerical simulations [ijj. Figure shows the dimensionless power density as a function of nt,/n s h- Taking into 




FIG. 7: Plot of the dimensionless power density as a function of the ratio of the bulk plasma density to the sheath density, 
taking into account (a) both Ei(x) and Eb - solid line; (b) only Eb - dashed line; and (c) no electric field inside the plasma - 
dotted line. 

account the electric field in the plasma (both Eb and Ex) reduces the total power deposited in the sheath region. 
Interestingly, taking into account only the uniform electric field Eb gives a result close to the case when both Eb and 
E\ are accounted for. The electric field E\ redistributes the power deposition from the energetic electrons to the low 
energy electrons, but does not change the total power deposition (compare lines (a) and (b) in Fig and Fig- El - 
Therefore, the total power deposition due to sheath heating can be calculated approximately from Eq. J33J), taking 
into account only the electric field E b . This gives 



-mV? / u 2 



u n b , . 

B(u - u sh ) - 1 

u n sh 



^du. (31) 
du 



The result of the self-consistent calculation of the power dissipation in Eq. I|31|) differs from the non-self-consistent 
estimate by the last term in Eq. i|31|) . which contributes corrections of order n s h/nb to the main term. 

A future development should provide a self-consistent analysis of a more realistic, nonuniform, and self-consistent 
ion density profile ni(x). Such study has been currently performed for inductively coupled discharges only. 



III. PENETRATION OF THE RF ELECTRIC FIELD INTO AN INDUCTIVELY-COUPLED PLASMA 

Low pressure inductively-coupled rf discharges are often operated in the non-propagating regime, when the driving 
rf field penetrates into plasma only within a skin layer of width S near the antenna, i.e., exhibits a skin effect. Not 
only the rf field, but, in this case, also the resulting induced electric current is concentrated near the surface of the 
plasma. Depending on the local, or non-local nature of the relation between the electric current j induced in plasma 
and the rf electric field E, the skin effect is called normal, if the dependence of the current on the electric field is 
local, or anomalous, if the dependence of the current on the electric field is nonlocal |29|. 
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To differentiate between the two regimes of the skin effect, it is convenient to introduce the nonlocality parameter 
[3l| A = (A/$o) 2 j where A = vx/{^ 2 + v 2 ) 1 / 2 is the effective electron mean free path and 

^° = 71~, 2 1 2M/4 ( 32 ) 

is the depth of the normal skin effect. The parameter A 



2 ? 

a = TLJP 

u 2 c 2 (i + v 2 /u> 2 y/ 2 y ' 

is a fundamental measure of plasma current non-locality. In the local limit A « 1, the effective mean free path is 
small compared with the skin depth A <C 5, and the current density at a particular point in space can be considered 
as a function of the electric field at the same point j(x) = <r(x)E(x) (Ohm's law). In the opposite limit A > 5, the 
mean free path exceeds the skin depth A> j, the relation between the current and the field j(x) = J ct(x, x')E(x')dx' 
is no longer local, because the conductivity cr(x, x') has a spatial dispersion. 

The penetration of the rf electric field into the plasma is described according to Faraday's and Ampere's laws 

VxE = -if, (34) 
c at 

„ 19D 41. 

VxB = -— + — j. 35 

c at c 

For a transverse harmonic wave in one-dimensional geometry E y (x)e~ lu)t , the Faraday's and Ampere's laws give 

d 2 uj 2 \ „ 4niu; 



dx 2 ' ^ 

where the current j is the plasma electron current j y = j ey (the ions are considered stationary), which has to be 
calculated making use of the electron kinetic equation, similar to the case of the penetration of the longitudinal wave 
into the plasma described in the previous section. 



A. Normal skin effect 



In the limit of the normal skin effect (A <C 1), the electron thermal motion can be neglected. The electron flow 
velocity V ey may be obtained from Newton's law taking into account the drag force due to the electron neutral 
collisions, 



d_ 

dt 

This gives for the electron current (j ey = —en e V ey ) the Ohm's law relationship 



m "Hl^ey = —eEy ~ vV ey - (37) 



where 



j e (a:) = a e E(x), (38) 



(39) 



m{v — ilu) 

The plasma current density is proportional to the electric field at the same point of space with a proportionality coeffi- 
cient that is the complex conductivity of the cold plasma. Substituting Ohm's law Eq. I|38|l with plasma conductivity 
from Eq. I|39|l into Eq. (|36f) gives the solution of the wave equation 

E y = E y0 e~ ax , (40) 

where a = \f— Amuja e / c 2 . Here, we neglected small terms associated with the displacement current in the limit 
u <C u) p , which is valid for the most plasma parameters in ICP discharges. The electric field can be equivalently 
expressed as 

E y (x, t) = E y0 e~ cos{e/2)x/s " cos[tjt - sm{e/2)x/S } 1 (41) 
where 5q is the normal skin depth in Ea. l|32l ). and e = arctan^/w). 
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B. Anomalous skin effect 



The case of anomalous skin effect (A > 1) for low-pressure inductively-coupled plasmas is more complicated com- 
paring to the case of normal skin effect, and requires a more elaborate mathematical and numerical treatment to 
uncover its intrinsic complexity. In the limit A ^> 1, the electron mean free path is large compared with the skin 
depth, and the electron current is determined not by the local rf electric field (Ohm's law), but rather is a function of 
the whole profile of the rf electric field over distances of order A. Therefore, a rather complicated nonlocal conductivity 
operator has to be determined for the calculation of the rf electric field penetration into the plasma. 

In the case of a uniform plasma, the Vlasov and Maxwell equations can be solved by applying a Fourier transform 
|28|. For a transverse harmonic wave in one-dimensional geometry E y (x)e~ lut , a spatial Fourier harmonic of the 
current j y kexp(~ikx) simplifies to become |2fll8f]j 

" l£k Z {lk) E >- (42) 

Details of the solution are given in Appendix C. The electric field profile is given by the inverse Fourier transform of 
Eq.® 

2? ui f°° e ikx 

Here, / is the surface current in the antenna and et{uj,k) is transverse plasma permittivity, which for a Maxwellian 
EEDF is given by 

Note the module sign as an argument of the plasma dispersion function. It reflects the proper symmetry of the 
continued electric field profile into semi-plane x < and also the proper pole position of the plasma dispersion 
function |2j, |9j]. Neglecting the module sign results in erroneous results. 

The solution for the electric field Eq.l^SJ has been described in many reviews and textbooks 0, E3, El> Ell • Here, we 
only focus on a property of the solution (|48|l not commonly acknowledged in the literature. 

In the limit A > 1 or i < w^/w, the plasma dielectric function can be substituted by its limiting value at small 
arguments Z ~ ii/ft. Introducing the anomalous skin depth 



c 



uj v \ujcJtt 



and substituting Z ~ iy/ir into Eq. (|44|l and into Eq. (|43|l gives 



2iu 



E v (x) = —I I 5 JTHJ^ dk. (46) 



The integral in Eq. (|46[) cannot be calculated analytically, but it can be transformed into an integral in the complex 
k plane by substituting \k\ = y/W. The contour of the integration should encompass branch point of the function 
y/W and has to come around the imaginary A:— axis. This gives 



(i V3 + l) ex (_xj2 \ + Eq 

371 V da J 3 7l ^ V &a 

UE 0n f 00 £exp(-xZ/6 a ) 



E y (x) = E ^i^ exp ( ) + ^ exp ( - f ) (47) 



7T7i Jo 1 - e 



d£. (48) 



where 71 = 2(v3 + i)/3v3, 72 = (1 — iv3)/2 and Eq is the electric field at the plasma boundary at x — 0, P stands 
for principal value of the integral. The last term represents the contribution of the integral around the imaginary 
k— axis and the exponential terms originate from the poles. The electric field at x = can be calculated analytically 

Aiul -k(\/Z + i)8 a 

E ° = —2 P72 ' ( 49 ) 
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From Maxwell's equations it follows that the magnetic field near the coil is -B|o+ = 2wl/c. Correspondingly, the 
derivative of the electric field at the plasma boundary is 



dE y , 2iritu 
d.r 

The characteristic decay length of the electric field can be introduced as |33j, |3 

E, 



Vi x =o=-^I- (50) 



~(l + i/Vt)6 a . (51) 



—dEy/dx 3 

The electric field profile from Eq. (|47|l is compared in Fig. [S] with the exponential profile 

E y {x) = E exp[-xRe(l/S s )] . (52) 
A more conventional plot of the amplitude and phase of the electric fields is shown in Fig. [5] 




x/S x/S 

a ; 



FIG. 8: Plot of the rf electric field as a function of the normalized coordinate x/S a - The solid curve corresponds to the solution 
in the limit A = vt(j-> p /clo = oo; dashed line - A = 93 (plasma parameters n = 10 11 cm -3 , T e = 3 eV, f = 1 MHz). The dotted 
and dash-dotted lines shows the skin approximation in Eqs. 151H and 1521 : (a) real, and (b) imaginary part of the electric field. 




FIG. 9: Plot of the rf electric field and electron current as a function of the normalized coordinate x/5 a - The same profiles as in 
Fig- HI shown are (a) amplitude, and (b) phase with respect to the phase of the electric field generated by the field in vacuum. 
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C. Spatially averaged electric field, J°° E y dx — > in the limit of a strong anomalous skin effect A — > oo. 

The most apparent difference between the anomalous skin effect and the normal skin effect is that the amplitude 
of the rf filed is non-monotonic in the limit of anomalous skin effect and monotonic (exponential) for the normal skin 
effect. Moreover, in the case of the extremely anomalous skin effect, in the limit A ^> 1, the spatially averaged rf 
electric field tends to zero 

/•OO 

Eydx -> 0, A -> oo. (53) 

In other words, the phase of the electric field changes by tt inside the skin layer, see Fig.[5Jb). The spatially averaged 
electric field is given by the Fourier component at k = 0, i.e., 



POO 

/ Eydx = nE(k = 0). 
Jo 



(54) 



Substituting the Fourier component of the electric field from Eq. (|43[1 into Eq. I|54|) gives 

f°° , . j 2iriujl 1 

/ Ey{x)dx = - — 

Jo c 1 ui 2 e(lu)/c 2 

and 



00 . . , 2niu;I 1 , , 

Ey{x)dx = — - (55) 



^E v {x)dx _ 33/M/3 

\E \6 a ~ - AJ73-' (56) 

From Eq. 1(56(1 it is evident that as the nonlocality parameter tends to infinity, the averaged electric field tends to 
zero. This property of the electric field profile is consistent with nonlocality of the electron current. The electric field 
profile and the current profile are coupled to each other by Eq. (|36|) . Therefore, the main part of the current and the 
electric field should decay on distances of order 5 a , see Fig. However, if the electric field profile has a non-zero 
average, the fast electrons will pick up a velocity kick from the skin layer and will transport the current over distances 
of order vtIuj 3> 5 a , where the electric field vanishes. This would contradict Maxwell's equations. Therefore, the 
zero average of the electric field is necessary and an important property of the electric field profile in the limit of the 
extreme anomalous skin effect A — > oo. 
The penetration length is defined in textbooks [33l l34| as 

X E = Io°° E ^ dx . (57) 

From the above discussion it follows that this definition is confusing, because in the limit of the anomalous skin effect 
the above defined penetration length is A# <C 5 a and is not a good measure of penetration length of the electric field. 
A better definition would be 

fn°° \E v (x)\dx , „ 

Vl= ° • (58) 

In the limit of the strong anomalous skin effect, i.e. A> 1, numerical calculation gives 

A| B | = 1.64£ a . (59) 

From Fig. |5| it is evident that in the region x < 28 a the amplitude of the electric field can be approximated by the 
exponential profile in Ea. (|52|l with the decay length 

5 < - Reim = I 6 - m 

Note that the penetration length defined by Eq. 1)58(1 . Ai£i is nearly twice as large as the initial decay length of the 
electric field amplitude near the plasma- wall boundary S e . This is due to the pronounced long tail in the profile of 
the electric field. 
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Similarly, if we introduce the penetration length of the current 



Jo°° \iv( x )\ dx 
liol 



% = J0 , ■ (61) 



numerical simulation gives 



A b1 = 1.87* « A| E |. (62) 



This result contradicts to claim of Refs. [331 l35fl. that the magnetic field and current penetration lengths are much 
longer than the electric field penetration length. This claim is the result of an inaccurate definition of the penetration 
length. 

In an attempt to reduce the phenomenon of the anomalous skin effect to the normal skin effect, many authors 
have substituted the correct profile of the electric field in Eq. I|47|l by an exponential profile Eq exp(— x/S e ) with some 
fitting procedure for S e |3a.l37ll3q . By doing so, the property of the electric field in the limit of anomalous skin effect 
in Ea. (|53|l is violated. This leads to overestimation of the electron heating 0. Under the conditions of the anomalous 
skin effect vt S> S a LO, electrons acquire a velocity kick 

2e f°° 

Av y = / E y (x)dx. (63) 



If E y (x) satisfies the condition in Eq. (|53|l . the electron velocity kick after passing through the skin layer is much 
smaller than in the case of an exponential electric field profile, which does not satisfy the property J °° E y dx — ► 0, as 
A — > oo. 



D. Analytical separation of the electric field profile into an exponential part and a far tail. 



Consider an exponential profile of the electric field in a plasma 

E y (x, t) = E y o exp (—kpX — iwt) , 
where k p is a real positive number. The velocity perturbation in this electric field becomes 



e ft 

Av y (x,t) = / d,TE y [x(T),T]. 



(64) 



(65) 



The velocity kick Av y can be separated into a purely exponential part and a non-exponential part. Substituting the 
electron trajectory x(r) = x — v x (t — r) for v x < gives 



Av y (x,t) 



Ey 



m —k p v x — ioj 



— exp (—kpX — iut) . 



(66) 



For v x > 0, the velocity acquired by an electron can be represented as the difference between the velocity kick acquired 
after a full pass through the skin layer and the contribution from the part of the skin layer [x; oo], i.e., 



Av y (x,t) 



dTE y [x(T),' 



(67) 



The second part of the integral (Awfj) in Eq. I|67|l gives an exponential profile for the velocity kick, similar to Eq. I)66|) 



Av y (x,t) 



E 



m —k p v x — ilu 
The first part of the integral (Av y n ) in Eq. I|67|l gives 



exp (— k p x — iut) , v x > 0. 



(68) 



Av y n = A^e-M*-*/^ 



(69) 
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Here, Ao™ is the velocity kick acquired during the pass through the entire skin layer. The time t — x/v x corresponds 
to the moment the electron collides with the wall. 

Substituting Av y (x,t) from Eqs. Ij66(l and l|68l ) gives for the exponential part of the current 

8f 

Av v —v y dv, (71) 

f y = -E y0 e-^ k ^ [ / . |^ % dv, (72) 
y m J —k p v x — iuu ovy 

or, after integration, the exponential profile of the current becomes 

fy = -E y0 e-^-^-^-Z f-^-X . (73) 
y m k p V T \k p V T J 

The asterisk denotes the complex conjugate. Note that Eq. H73[) can be derived from Eq. (|42[) with the substitution 
k = ik p and by accounting for the following property of the dispersion function 

Z(C) = ~Z(-0*. (74) 

The exponential part of the profile should satisfy Maxwell's equation (|36|) . This gives an expression for k p 

Note that because Z in Eq. I|75|l has only purely imaginary and positive parts, k p is a real positive number, as it was 
assumed to be. 

The non-exponential part of the electron velocity kick in Eq. I|69|) generates a non-exponential part of the current 
profile, which decays over a spatial scale of order Vr/ui due to the phase mixing, as the phase of the velocity kick 
u(t — x/v x ) in Eq. (|69J) is different for electrons with different v x . The current and electric field profiles are essentially 
non-exponential, similar to Eq.Q) for longitudinal velocity kicks, as discussed above. 

Details of the exact analytical calculation of the electric field profile separation is given in Appendix C. Applying 
a procedure similar to that of Landau's treatment 2] for the longitudinal electric field, the integral in k— space in 
Eq. 1431) can be separated into an integral over an analytic function in the region k € [— oo, oo] and an integral over 
some non-analytic function in the region k € [0, oo]. To do so, the plasma permittivity has to be analytically continued 
from the real axis k < 0, Irak = 0, into the complex k -plane, see Appendix C for details. The first integral can be 
readily calculated using the theory of residues. In the upper half-plane of the complex fc, there exists only one pole of 
the analytically continued function of the plasma permittivity continued from k < 0. The value of the pole is equal 
to ikp, given by Eq. (|75|l . 

In the limit u> ^ k p Vr, Z(Q) — — l/£, where £ = iui/kpVr- Substituting this value for the plasma dielectric function 
into Eq. I|75|l yields k p = t^p/c, i.e., the normal skin layer length l/k p = So in Eq. H32|) for v <C uj and lu <C lo p . Figure 
1 101 shows the profile of the electric field for the same typical ICP parameters: plasma density n = 10 11 cm -3 , electron 
temperature T e = 3 eV, and discharge frequency / = 13.56 MHz. Shown are the exact electric field profile E y (x) 
calculated according to Ea. H43|) . the exponential part of the electric field 

E V p{x) = E exjp(-kpx) (76) 

with k p from Eq.JTSJ), and the difference of the two 

E yt (x) = E y (x) - E yp (x), (77) 



and the asymptotic calculation for E yt {x) in Eq.jC22|) E yst d(x). For these plasma parameters the skin effect is neither 
normal nor anomalous: ui/kpVr = 1.52. Notwithstanding the fact that the parameter ui/kpVr is of order unity, the 
main part of the electric field is close to the exponential profile in Eq. (|76(l with k p from Eq. I|75|l . E y (x) ~ E yp {x). 
As evident from Fig. I1UI the non-exponential part is small, E yt (x) <ti E yp (x), everywhere where the electric field is 
substantial, or up to distances five times of skin depth, for x < 5/k p — 7.5Vr/w. The tail of the electric field profile 
for x > 7Vt/uj is non-exponential and dominated by E yt (x). 

In the limit of the anomalous skin effect u>/k p Vr -C 1, Z(Q) — iy/ir, where £ = iui/k p Vr- Substituting this value for 
the plasma dielectric function into Eq. I|75|) yields k p = 1/S a , which is very close to the skin impedance approximation 
in Eq. I|52|l which corresponds to k p = 9/85 a -a 12 % difference. As a result, the exponential profile in Eq. I|76|) 
approximates well the exact profile of the electric field over distances within a few skin depths even in the limit of the 
strong anomalous skin effect, as is evident in Fig. 1111 However, the non-exponential part E yt (x) dominates E yp {x) at 
x > Vt/lo in accord with the requirement in Eq. I|53|) . 
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FIG. 10: Plot of the rf electric field as a function of the normalized coordinate xuj/Vt for plasma parameters n = 10 11 cm -3 , T e = 
3eV, / = 13.56 MHz. Shown are (a) the amplitude and (b) the phase. Solid lines show the exact electric field profile E(x) 
calculated according to Eq. 1431 : dashed (red) line, the exponential part of the electric field E p (x) = Eq exp(—k p x) with k p from 
Eq. I|75|l : dotted line (green), the difference of the two Et(x); and, chain (cyan) line, E s td(x) shows the asymptotic calculation 
for Et in Eq. IC22I . Subscript y is ommitted in the electric field. 




FIG. 11: Plot of the rf electric field as a function of the normalized coordinate xui/Vt for plasma parameters n = 10 n cm -3 , T e = 
3 eV, / = 1 MHz, similar to Figs. |SJ and |S] Shown are (a) amplitude and (b) phase. Solid lines show the exact electric field 
profile E(x) calculated according to Eq. H43H : dashed (red) line, the exponential part of the electric field E p (x) = Eo exp(—k p x) 
with k p from Eq. 117511 : dotted line (green) the difference of the two Et(x); chain (blue) line represents the limiting case of strong 
anomalous skin effect A — > oo E apP r(x), and dashed and double dotted (chain) line shows E st d(x), the asymptotic calculation 
for E t in Eq. 1C221 . Subscript y for the electric fields is omitted. 



E. Surface impedance 



An important plasma characteristic is the surface impedance, which is given by the ratio of the electric field to the 
rf magnetic field or the coil current at the plasma boundary 3] 



where 



^=§l*=o, (78) 



B|*=o = ^-1 (79) 
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is the magnetic field near the antenna. The total power P deposited per unit area into the plasma is determined by 
the energy flux dissipated into the plasma or the time-averaged Poynting vector 

P=<S x >=~Re(EB*). (80) 

Substituting the electric field from Eq. l(78|l and the magnetic field Eq. (|79|) into Eq. ||%U|> relates the power to the real 
part of the surface impedance 

P = ^-I 2 ReZ. (81) 
2c 

The imaginary part of the surface impedance describes the plasma inductance. 

The surface impedance can also be used to estimate the penetration length in the surface impedance approximation 
given by Eq. (51) • Substituting the electric field from Eq. 17811 and the magnetic field Eq. (|79|) into Eci. (|51J) relates 
the penetration depth and the surface impedance 

r7 

6 S = —. (82) 
ioj 

The surface impedance can be calculated making use of Eq. j|, i.e., 

/OO ^ 
k*-u*e t (o,,k)/<? dk ' (83) 

which requires numerical integration. On the other hand, we can use the results of the previous subsection that the 
main part of the electric field is an exponential function in Eq. I|76|) with k p given by Ea. (|75() . From Eq. I|82|) . the 
imaginary part of the surface impedance can be obtained substituting S s — l/k p 



Z„ = » . ,84, 



A pure exponential profile yields only the imaginary part of the surface impedance. The real part of the impedance 
can be calculated by computing the power dissipated by electrons from the skin layer |36j ] 

P=^|^/|A^| 2 dv, (85) 

where Ao™ is the velocity kick acquired by an electron after passing through the skin layer, which is given by Eq. (7UJ|. 
Here, m(An™) 2 /4 is the temporal average of the electron energy change in the skin layer and v x f is the electron flux 
on the wall. Equation (|81|l becomes 

2 o , „ , o f / knVr,. 



ReZ p = ^ J dVx . (86) 

Because the imaginary part of impedance is large compared with its real part, only the imaginary part can be included 
on the right hand side in Eq. i|86|) . Figure IT2"1 shows the real and imaginary parts of the surface impedance versus the 
discharge frequency calculated exactly, i.e., making use of Eq. l|83|l . and approximately from Eqs. (|84|) and l|86l) . Also 
shown at the top of this figure, is the ratio of the actual skin depth 5 = l/k p from Eq. I|75|) to the normal skin depth 
calculated in the cold plasma approximation Sq given by Ea. (|32l) . From Fig. ll2l it is evident that within 50% accuracy, 
the impedance calculation can be based on the exponential profile in Eq. (|76|) for discharge frequencies higher than 1 
MHz 37] . However, for lower frequencies, the assumption of purely exponential profile leads to overestimation of the 
electron heating and plasma resistivity up to a factor of 3 for / <C 1 MHz, see Fig. 1121 This is because the important 
property of the electric field profile under the conditions of strong anomalous skin effect in Eq. I|53|l is being violated. 
Note that at these low frequencies taking into account a small but finite collision frequency or nonlinear effects may 
be important. 



F. Anomalous skin effect for an anisotropic electron velocity distribution 

The anomalous skin effect in a plasma with a highly anisotropic electron velocity distribution function (EVDF) 
is very different from the skin effect in a plasma with the isotropic EVDF. In Ref. |40| an analytical solution was 
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FIG. 12: Plot of the real and imaginary parts of the surface impedance versus discharge frequency calculated exactly making 
use of Eq. 1831 and approximately using Eqs. 1841 and 1861 in the limit of collisionless plasma v <C oj. Also shown is the ratio 
of the actual skin depth S = l/k v given by Eq. 1751 to the skin depth calculated in the cold plasma approximation So Eq. 1321 . 
(top). 



obtained for the electric field penetrating into plasma with the EVDF described by a Maxwellian with two temperatures 
T y T x , where y is the direction along the plasma boundary and x is the direction perpendicular to the plasma 
boundary. Under the conditions 



VTy C 

— i > — ; oj p > u, 
OJ aj v 



(87) 



the skin layer was found to consist of two distinct regions of width of order vt-x/u and vt v /u, where vx x ,y = \/T x ,y/m 
are the thermal electron velocities in x and y directions, and lo is the incident wave frequency. The calculation is 
based on Ea. (|43|l . where the dielectric permittivity has to be modified for an anisotropic EEDF to become 



e t (uj, k) 



1 



1 



1 



-Z 



(88) 



V2v Tx k 

In the case of anisotropic EEDF under conditions in Eq. (|87() , the integral in Eq. (|43(l has two poles and the integration 
over the branch point k = does not contribute. As a result, the profile of the electric field is a sum of the two 
complex exponents: 



E(x) ~ —B(0) 

L0 r , 



UpVTy 



exp(ik pl x) + ^T x /T v exp(ik P2 x) 



where k p \ is given by 



and k P 2 is given by 



CO 



VTy 



k P 2 — —\/T y /T x 



(89) 



(90) 



(91) 



c v J 2V2v Tx ' 

The profile of the electric field is shown in Fig. ^| The skin layer contains multiple oscillations of the electric field, 
in striking contrast to the case of isotropic EEDF. 
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FIG. 13: The electric field in the plasma with vt v = 0.1c, u) — 0.01oj p , T y /T x = 50. The solid line shows the real part 
of the electric field profile obtained from the full solution. The dashed line corresponds to the smooth part of the solution 

~ exp(— LOx/VTy)- 



IV. CONCLUSIONS 



We showed that electrons can transport the plasma current away from the skin layer due to their thermal motion 
over distances of order Ut/w. As a result, the width of the skin layer increases when electron temperature effects 
are taken into account. The anomalous penetration of the rf electric field occurs not only for the wave transversely 
propagating to the plasma boundary (inductively coupled plasmas), but also for the wave propagating along the plasma 
boundary (capacitively coupled plasmas). It was shown that separating the electric field profile into exponential 
and nonexponential parts yields an efficient qualitative and quantitative description of the anomalous skin effect. 
Accounting for the non-exponential part of the profile is important for the calculation of the electron heating and the 
plasma resistivity. For example, the assumption of purely exponential profile leads to overestimation of up to a factor 
of 3 in the electron heating for / <C 1 MHz, see Fig. 1121 

Here, we considered only plasmas with a Maxwellian electron energy distribution function. However, in low pressure 
rf discharges, the EEDF is non-Maxwellian for plasma densities typically lower than 10 10 cm -3 |32|. The nonlocal 
conductivity, and plasma density profiles and EEDF are all nonlinear and nonlocally coupled [3jJ . Hence, for accurate 
calculation of the discharge characteristics at low pressures, the EEDF needs to be computed self-consistently [4lll42l 
The effects of a nonMaxwellian EEDF, nonlinear phenomena, the effects of plasma non- uniformity and finite 
size, as well as influence of the external magnetic field on the anomalous skin effect will be reported in the second part 
of the review [4o| . 
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APPENDIX A: ANALYTICAL DERIVATION OF THE CURRENT PROFILE DRIVEN BY VELOCITY 

KICKS NEAR THE PLASMA BOUNDARY 



Consider that electrons acquire a velocity kick near the boundary, in the direction perpendicular to the boundary 

dv x = AUcos(wt). (Al) 

The electron velocity at a distance x from the boundary will be determined by the exact moment of the collision with 
the boundary at a time t — x/v x . The electron current in the plasma is given by integration over all electrons with a 
distribution function f(v x ) 

/>oo 

j(x,t) — eAV / f(v x ) cos(wt — ujx/v x )dv x . (A2) 
Jo 

For a Maxwellian distribution function f(v x ) = nQe~ v *l VT /vt^/k the current in Eq. I|A2(1 takes the form j(£,t) = 
joA(£) cos[wt — </>(£)] , where j'q = enoAV and A and <f> are the amplitude and phase of the current, respectively, and 
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£ = lux/vt- The functions A and are shown in Fig. 1. In the limit £ 3> 1, the integration in Eq. i|A2(l can be 
performed analytically making use of the method of steepest descend [J 



Jo 



Re e 



(A3) 



where and s = v x /vt- The integral in Eg. (|A3|> can be calculated in the complex s plane. The stationary phase point is 
given by d(— s 2 + i£ / s) / ds = or s 3 = — ££/2. This gives the stationary point s = (— i^/2) 1 ^ 3 . In the neighborhood of 
this point, the function in the exponent can be expanded as a Taylor series, — s 2 + i£/s = — Sq + i^/sq — 6(s — so) 2 /2 = 
— 3sq — 3(s — so) 2 . Integration of the Gaussian gives / °° e - 3 ( s ~ s o) d s = y^n~/3. Substituting this into the integration 
m Eq. JSHJl yields: 



t) = ^=Re (exp [-iojt - 3 (-^C/2) 273 ] ) . 



Substituting (— i) 



2/3 



(e-*r/2) 



2/3 



e - l7r /3 _ cos (^-/3) _ i s in(7r/3) = 1/2 - V3i/2 into Eq. (jSiJ) gives 

JO __ ( , 2 /3 /o5/3\ CQS _ 3 V3 e 2/3 /2 5/3> 



,(£./) =^e, P ( -3£ 2 / 3 /2 5 / :! 



(A4) 



(A5) 



APPENDIX B: ANALYTICAL DERIVATION OF THE LONGITUDINAL RF ELECTRIC FIELD 
PROFILE NEAR THE PLASMA BOUNDARY (E || k) 

The analytical solution for a longitudinal rf electric field involves solving the Vlasov equation for the electron 
velocity distribution function (EVDF) F 

dF dF e _ OF n 

-K7+Vxj> E x —=0, Bl 

ot ox m ov x 



^ = Aire [ n t - I Fdv x ) . (B2) 



together with the Poisson equation 

dx 

In the linear approximation, the EVDF can be split into two parts 

F(t, x,v x ) = f (v x ) + f(t,x,v x ), (B3) 

where fo(v x ) describes EVDF of a uniform plasma with uniform ion density n e = tii = uq and f(t,x,v x ) is EVDF 
due a wave perturbation. Substituting Eq. I|B3I) into Eqs. <|B1|) and (|B2|) yields the linearized Vlasov-Poisson system 
of equations 

-wr + v xa E x — = -vf, (B4) 

at ox m dv x 



dF f°° 

= -4ne / f(v x )dv x . (B5) 
dx J_ QO 

In the first equation l|B4|) . the small collisional term with the collision frequency v <C oj is taken into account. In 
Ref. [2| Landau solved the linearized Vlasov-Poisson system making use of the Laplace transform for a semi-infinite 
plasma x > 0. However, it is more convenient to apply a Fourier transform to an infinite plasma by artificially 
continuing the EVDF and the electric field in the semi-plane x < . Electrons moving with v x < reflect from the 
boundary x — and change their velocity to — v x . This gives the boundary condition for the Vlasov equation in the 
semi-plane x > 



f{t,0,v*) = f(t,0,-v x ). 



(B6) 
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Instead of considering problem in the semi-plane x > with the boundary condition in Ea. (|B6|) . we can consider the 
entire plane x € [—00, 00] by artificially continuing the electric field into the semi-plane x < 0. The Vlasov equation 
is symmetric with respect to a change in variables according to the substitution 

v x — * —v x , x — > -x, E — > (B7) 

Therefore, electrons at x = with > 0, which are reflected from the wall can be represented as electrons which 
came from the semi-plane x < and interacted with the electric field 

E x (x < 0) = -E x (x > 0). (B8) 

As a result, the electric field has to be continued anti-symmetrically into the semi-plane x < 0. 

Now we can apply the Fourier transform for the Vlasov-Poisson system of Eqs. I|B4|) and i|B5(l . This gives for the 
components of the EVDF f k e~ iuJt+kx and the electric field E k e~ lut+kx 

-i(u + iv-v x k)f k --E k J^-=0, (B9) 
m dv x 



ikE k + 2E = Awe 



fkdv x 



(B10) 



Note that due the fact that the electric field is a discontinuous function, the Fourier transform of the derivative of 
the electric field dE/dx is ikE + 2Eg, where Eq = E(0) is the electric field at the right side (x > 0) of the plasma 
boundary. Substituting f k from Ea. (|B9() into (|B10|I yields 



E k = 

where eh(w, k) is the longitudinal plasma permittivity 

£{l (u J ,k) = l + ^- 
" nok 

Substituting a Maxwellian EEDF 



2E n 



1 



ik e\\(u>,k) 



1 



dfo 



fo = 



no 



' ttvt 



uj + iv — v x k dv x 



■exp(-w 2 /w|.), 



dv x . 



where vt = \/2T/m, into Ea. l|B12JI and after some algebra 0], we obtain 

f 



£||(w,fc) ~ 1 + 



/TTVT 



50 exp(-v 2 /vl) 
x v x k - u> - iv 



The last term on the right hand side can be expressed in terms of the plasma dispersion function 



HO 



1 



exp(— t 2 ) 



dt, Im(C) > 0. 



(Bll) 



(B12) 



(B13) 



(B14) 



(B15) 



The dispersion function Z(Q in the form of Eq. l|B15p is only defined for Im(Q > and is defined as an analytical 
continuation for Im(() < 0. For k > 0, in the limit v — > 0, 



'TTVT 



exp(-i; 2 /i;|,) 1 
v x k — u> — iv x kVr/tjJ 



Z(aj/kV T 



(B16) 



For k < 0, the imaginary part of the (lo + iv)/k is negative and we have to transform the integral i|B14fl so that the 
pole v xp = (lu + iv) I k lies in the upper plane of the complex velocity. This can be achieved by substitution — v x — > v x 
, which gives for k < 



exp(-w 2 /«T) , 
-dv x 



1 



TTVT J -00 v x \k\ - u - iv 



\k\vT/>JJ 



Z(u/\k\vr). 



(B17) 
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As a result, 



en(uj,k) ~ 1 + 



2w 2 
k 2 v\ 



\kV T /uJ 



■Z{\u/kV T \) 



(B18) 



Note that because the function fg is symmetric with respect to the substitution —v x — > v x , s(lo, k) is symmetric with 
respect to the substitution — k — > fc. Correspondingly the symmetry of the electric field in Ea. i|B8|l is preserved. 
The electric field profile is given by the inverse Fourier transform of Eq. ( IB11|) 



Z7T 



2E e ikx 
ik £u(u),k) 



dk. 



(B19) 



In the limit x — > oo, E(x) — > Eq/s, where e = s(u>,0). This is in accord with the conservation of the total current 
in the one-dimensional geometry. The total current is the sum of the displacement current and the electron current, 



1 dE x 



+ Je=I(t). 



(B20) 



The total current conservation follows from the combination of the Poisson equation and the charge continuity 
equation. Indeed, taking the time derivative of the Poisson equation and making use of the charge continuity equation 
gives 



d 

v-— /•:, • i-V-./. 
ot 



0. 



(B21) 



In one-dimensional geometry it can be integrated with a constant of space - the total current carrying through the 
plasma I(t), which gives Ea. i|B20|l . For a harmonic electric field considered here, Ea. lB20|l gives 



iluE t — Attiuj 



e- 1 

47T 



E x 



-iuiEo. 



(B22) 



Here, we account for the relationship between the plasma conductivity (J e — aE) and the plasma dielectric function 
e = 1 + 4na/(-iu>). Ea. (|B22|) gives 



E x (x — > oo) = E /e. 



(B23) 



The same result can be obtained from Ea. (|B19|l after substituting e(u>, k) — > e(cu, 0) and integrating. Thus, the electric 
field in the transition region is given by 



1 /°° 2B 
E x (x)-E /s=- 



1 



1 



ik \su(uj,k) e\\(u>,0) 



3 ikx dk. 



(B24) 



The dielectric function in the form given by Eq. (|B14 |) is not an analytic function of k. To apply the theory of residues, 
Landau proposed to split integral into two parts |2| according to 



ik \E\{uj,k) e(o;,0) 
2E 



2ir Jo ik \Ei(u>, k) e\\(u>,k) 



3 lkX dk 



e lkx dk, 



(B25) 



where 



e x (u,k) = 1 + 



LU 2 k 2 



1 - 



kvT /to 



Z( — Lu/kVT) 



(B26) 



The first integral can be calculated by moving the path of integration into the complex A:— plane and applying the 
theory of residues. For to <C uj p , e < and there is only one pole ei(u>, k) = in the upper half-plane |2|. It corresponds 
to the usual screening with the Debye length. In the limit k ~ lj p /vt, Z(\u>/kvT\) ~ 1 and Z {\uj / kvrY) / \kvT / o->\ <C 1, 
which gives 



£\(UJ, k) ~ 1 + 



H 

k 2 v\ ' 



(B27) 
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Calculation of the first term in Ea. i|B26|l gives Eq exp(— x/a), where a is the Debye length a = vt/V^^p- Therefore, 

(B28) 



E x {x) = E /e + E exp(-x/a) + — / — - — — 

27r J Q ik \ eiiuj, k) e\\{u, k) 



e lkx dk 



For Im(k) = 0, Z(-oj/kv T ) = -Z(uj/kv T )* B and 

El(w,fc) = l + " 



v\k 2 



kvT /u> 

Substituting Eq. (|B29fl into Eq. [|B28fl gives for the last term E t (x) 



where Q 



Z(uj/kvT) 



= 4Eo^4 /- 1 Im[Z( U /kvT)) eikxdk 
ir uf, y fc 4 ei(w, fc) £||(w,fc) 



Jm[Z(C)] = ^exp(-C 2 ). 



(B29) 



(B30) 



(B31) 



The last integral can be calculated analytically only in the limit x ^> vt/u by applying the method of steepest 
descend. In this limit, k -C Vt/oj, £i(uj, k) s» e(u>, k) » e and 



/ -|r exp(ifcr - uj 2 /k 2 v^)dk ~ ^^(xA w ) 2/3 A 1 . 
Jo v 3 



exp 



2/3 

C ( f~ ) _ i7F / 3 



(B32) 



where c = 3(— 1 + i\/3)/4, and A w = vt/V^cj is the phase-mixing scale. 

Substituting Eq. l|B31|) into Ea. HB30|) and making use of Eq. <jB32|) yields at x » A w 



2E ^ f x \ 2/3 

w w vw UJ exp 



A. 



2/3 



17T/3 



(B33) 



The plots of amplitude and phase of the electric field profile E t (x) given by Eq.{B30j and the approximate analytical 
result Eq. (|B33|) are shown in Fig. [5] 

APPENDIX C: ANALYTICAL DERIVATION OF THE TRANSVERSE RF ELECTRIC FIELD PROFILE 

NEAR THE PLASMA BOUNDARY (E _L k) 



The analytical solution involves solving the Vlasov equation for the electron velocity distribution function (EVDF) 



F 



OF dF e 
dt dx m 



(E y +v x x B z ) 



dF 

dv v 



This equation has to be solved together with the Maxwell's equation yielding 



dx 2 c 2 v 



-tt E v = 



[J + IS(x)} , 



(CI) 



(C2) 



where I is the surface current. The plasma density is not perturbed in the transverse wave; therefore there is no need 
to solve the Poisson equation. In the linear approximation, the EVDF can be split into two parts 



F(t,x,v) = f (v) + f(t,x,v), 



(C3) 



where fo(v) describes EVDF of an isotropic, uniform plasma with uniform ion density n e = Uj = no and f(t,x,v) is 
the EVDF due a wave perturbation. Substituting Eq. (|C3|) into Eqs. (|C1|) yields the linearized Vlasov equation 



dt x dx m v dv y 



(C4) 



25 



In Eq. I|C4|) . the small collisional term with collision frequency v <C uj is taken into account. Similarly to the case of 
the longitudinal wave, we can consider the entire plane x € [—00,00] by artificially continuing the electric field in the 
semi-plane x < 0. The Vlasov equation is symmetric relative to the change in variables according to the substitution 

v x ~> ~v x , x -> -x, E y -> Ey. (C5) 

Therefore, electrons at x — with v x > which are reflected from the wall can be represented as electrons which 
came from the semi-plane x < and interacted with the electric field 

E y {x < 0) = E y (x > 0). (C6) 

As a result, the electric field has to be continued symmetrically into the semi-plane x < 0. 

Now we can apply the Fourier transform for Eqs. I|C4|I and (|C2|I . This gives for components of the EVDF f ke -^ t +' lkx 
and the electric field Ey k e~ iu ' t+ikx 

- i(uj + iv - v x k)fk - — E yk ^p- = 0, (C7) 

TO OVy 

^ + ^)E v ^^( Jk+ T } . (C8) 
Substituting f k from Eq. (|C7|) into (|C8|) with the current j k — —ej f k v y dw yields 

c 1 k 2 ~^e t {uj,k) 

where et(u), k) is the transverse plasma permittivity 



e t (u,k) = l + ^- r T7^ dv *- ( C1 °) 

no<jJ J -00 w + w - v xk ov v 



Substituting a Maxwellian EEDF gives Q 



£t{ ^ k) = 1 + ^»z(-?--). (Cll) 
oj z VT\k\ \VT\k\ ) 

Note that because the function /o is symmetric relative to the substitution — v x — > v x , s(co,k) is symmetric relative 
to the substitution —k^>k. Correspondingly, the symmetry of the electric field in Eq. (|C6|) is preserved. 
The electric field profile is given by the inverse Fourier transform of Eq. I|C9|) 

^»°i;L7-' f-^,^ (C12) 



Similar to the analysis of the longitudinal wave, we split the integral in Eq. I|C9|I into two parts 

E y {x) = E yp {x) + E yt {x), (C13) 

where 

2iriu> f°° p ikx 
E yv {x) = —l\ ,„ 2 „ a , t ,.^, 2 dk, (C14) 



k 2 — Lu 2 Eti(u>, k)/c 2 



and 



E (x) r [s t ^k)-e tl (u,k)]e^ 

£jyt{X > ~ c 2 c 2 J [k 2 - tu 2 e tl {u, k)/c 2 } [k 2 - u 2 e t (u, k)/c 2 ] ' 1 ' 

en {^k) = l-^^-z(-^-\. (C16) 
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Note that e t i(ui,k) = e t (ui,k) for fc < 0. 

The first part E yp (x) of the electric field can be calculated by evaluating the integral in the complex fc— plane. A 
pole of E yp (x)- k p i lies on the imaginary axis of the fc— plane. The dielectric permittivity is real and negative on 
imaginary axis of the fc— plane 



£ t i(w, kpi) = 1 1 



-F 



( 10 



ui 2 vrkp \VTkpJ ' 

where F(() = ImZ (i^) = y/n exp(y 2 )erfc(y) There is always a real value of k p as the root of 

kp = -uj 2 e a {uj, ik p )/c 2 . 
Applying the theory of residues, the integral for E yp {x) gives 



where 



2lLU 



2ttH 



c 2 2k p i — deti{uj,k) / dk lo 2 / c 2 



(C17) 



(C18) 



(C19) 



(C20) 



In the limit x S, the last term E yt {x) can be calculated making use of the method of steepest descend. Substituting 
fc 2 — u 2 Eti(uj, k)/c 2 in the denominator of the expression for E yt (x) by its limit = uj 2 /c 2 at fc — > 0, gives 



which yields 



E yt (x) 



E yst d(x) 



UJ 2 V T 



4j 2 7T V2 / x 
^JpVT \^ 



■ exp 



-1/3 



(—)' 

\vxk J 



exp 



+ ikx 



2/3 



dk 



c ( — ) - in/2 



(C21) 



(C22) 



where c = 3(— 1 + i\/3)/4 and A w = vt/V2^ 
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